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Abstract — The question of controllability is investigated for 
a quantum control system in which the Hamiltonian operator 
components carry explicit time dependence which is not under 
the control of an external agent. We consider the general 
situation in which the state moves in an infinite-dimensional 
Hilbert space, a drift term is present, and the operators driving 
the state evolution may be unbounded. However, considerations 
are restricted by the assumption that there exists an analytic 
domain, dense in the state space, on which solutions of the 
controlled Schrodinger equation may be expressed globally in 
exponential form. The issue of controllability then naturally 
focuses on the ability to steer the quantum state on a finite- 
dimensional submanifold of the unit sphere in Hilbert space - 
and thus on analytic controllability. A relatively straightforward 
strategy allows the extension of Lie-algebraic conditions for 
strong analytic controllability derived earlier for the simpler, 
time-independent system in which the drift Hamiltonian and 
the interaction Hamiltonia have no intrinsic time dependence. 
Enlarging the state space by one dimension corresponding to the 
time variable, we construct an augmented control system that can 
be treated as time-independent. Methods developed by Kunita 
can then be implemented to establish controllability conditions 
for the one-dimension-reduced system defined by the original 
time-dependent Schrodinger control problem. The applicability 
of the resulting theorem is illustrated with selected examples. 

I. Introduction 

Over the last two decades, quantum control has played an 
important part in theoretical and experimental progress toward 
the realization of laser control of chemical reactions and the 
development of quantum computers [1], [2], [3], [4], [5], [6], 
[7], [8], [9], [10], [11], [12], [13]. Essential to this contribution 
has been the integration of concepts and mathematical results 
from control engineering with the fundamental principles of 
quantum theory. 

Geometric control, a treatment of differential equations 
rooted in differential geometry, unitary groups, and Lie al- 
gebras, provides a natural mathematical basis for quantum 
control theory. Explicitly or implicitly, its elements [14] per- 
vade the manipulation of quantum states in both traditional 
and novel technologies. Indeed, the field of nuclear magnetic 
resonance (NMR) is largely concerned with geometric control 
of collections of interacting nuclear spins [12], [15], [16], [17]. 
Geometric control is also a key ingredient in the theory of 



quantum computation, figuring prominently in the works of 
Lloyd [18], Deutsch [19], and Akulin [20]. 

In particular, Lloyd [18] was among the first to establish that 
almost all quantum logic gates are universal. More precisely, 
if one has available a gate that can operate on two qubits, 
plus a single-qubit operation, then an arbitrary unitary trans- 
formation on the variables of the system can be performed 
with arbitrary precision by implementing a finite sequence of 
local operations. Clark [21] and Ramakrishna and Rabitz [22], 
[23] called attention to the close relationship between open- 
loop geometric quantum control methods and the application 
of quantum logic gates [19], [18]. 

Following Ref. [23], let us consider differential system 

-^l=AX(t)+Y,B t X(t) Ul (t), X(0) = I, (1) 

2—1 

which arises both in quantum computing and molecular con- 
trol. Here, X is a N x N unitary matrix (/ being the 
corresponding identity matrix), the matrices A and Bi, i = 
1, . . . , m are N x N skew-Hermitian, and the functions Ui(t) 
are controls. This equation is the law of motion of the 
evolution operators which govern time development of the TV- 
dimensional vector representing a pure state of the system in 
its iV-dimensional Hilbert space. A necessary and sufficient 
condition for Q to be controllable is that the set of all matrices 
generated by A, i = 1, . . . , m, and their commutators (i.e., 
the Lie algebra generated by A and Bi) equals the set of 
all N x N skew-Hermitian matrices. Additionally, when this 
condition is met, any X can be attained through some choice 
among the controls Ui(t) restricted to piecewise constant 
functions of time. In fact, the formulation adopted by Lloyd 
[18] in his universality proof corresponds to the special case 
A = and m = 2 of system Q. Already in the 1970s, 
Sussmann and Jurdjevic [24], [25] applied Lie-group theory to 
obtain rigorous results on controllability for finite-dimensional 
control problems corresponding to 0. 

Quantum computation has mostly concerned itself with 
the manipulation of discrete systems with finite-dimensional 
state spaces. However, the fundamental quantum observables 
representing position and momentum, and functions thereof, 
are continuous in nature. In view of recent developments 



in quantum error correction [26], [27], [28] and quantum 
teleportation [29], [30] of continuous variables, the potential 
of quantum computation over continuous variables warrants 
serious investigation, thus reopening issues of controllability 
on infinite-dimensional Hilbert spaces. Continuous quantum 
computers may in fact be able to perform some tasks more 
efficiently than their discrete counterparts. 

As early as 1983, Huang, Tarn, and Clark (HTC) [5], [31] 
proved a basic theorem on strong analytic controllability of 
quantum systems. This theorem explicitly embraces the case of 
quantum systems whose observables are continuous quantum 
variables acting on an infinite dimensional state space, but 
the essential finite-dimensional results may be extracted as 
special cases. Because of the difficulties caused by infinite- 
dimensionality and the unboundedness of operators, an ana- 
lytic domain in the sense of Nelson [32] was introduced to deal 
with domain problems [5], [31] and maintain key features of 
the application of Lie algebraic methods to finite-dimensional 
problems. 

Infinite-dimensional control systems have been widely if not 
systematically studied outside the quantum context. Brock- 
ett [14] addressed the problem of realization of infinite- 
dimensional bilinear systems. Sakawa [33] introduced a 
method for design of finite-dimensional controllers for 
diffusion systems with bounded input and output operators 
by using residual model filters. Keulen [34] designed infinite- 
dimensional Hoo controllers for infinite-dimensional systems 
with bounded input and output operators by using the solutions 
to two kinds of Riccati equations in an infinite-dimensional 
space. Based on gap topology, Morris [35] constructed finite- 
dimensional Hoo controllers for infinite-dimensional systems 
with bounded input and output operators. Morris [36] also 
showed that approximations of Galerkin type can be used to 
design controllers for an infinite-dimensional system. Costa 
and Kubrusly [37] derived necessary and sufficient conditions 
for existence of a state feedback controller that stabilizes a 
discrete-time infinite-dimensional stochastic bilinear system 
and ensures that the influence of the additive disturbance on 
the output is smaller than some prescribed bound. In Ref. [38], 
optimizability and estimatability for infinite-dimensional linear 
systems are investigated; also, a theorem on the equivalence of 
input-output stability and exponential stability of well-posed 
infinite-dimensional linear systems is established. In Ref. [39], 
the Hilbert-space generalization of the circle criterion is used 
for finite-dimensional controller design of unstable infinite- 
dimensional systems. There is also literature on absolute 
stability problems and open-loop stability problems in infinite- 
dimensional systems [40], [41], [42], [43], [44]. In addition, 
the spectral factorization problem plays a central role in 
designing feedback control for the linear quadratic optimal 
control problem in infinite-dimensional state-space systems 
[45], [46], [47], [48]. In contrast to this body of work, very 
little has been published on controllability for time-dependent 
infinite-dimensional quantum control systems. 

In the microscopic world ruled by quantum mechanics, 
most interesting phenomena involve change, and all real-world 



quantum systems are influenced to a greater or lesser extent by 
interactions with their environments. The environment changes 
with time, so the Hamiltonians used to describe these open 
quantum systems are explicitly time-dependent, as in Ref. [50], 
[49]. Tailored time-dependent perturbations are used to im- 
prove system performance [49] in high-resolution NMR spec- 
troscopy, where versatile decoupling techniques are available 
to manipulate the overall spin Hamiltonian [16]. Colegrave and 
Abdalla studied quantum systems with a time-dependent mass 
to investigate the field intensities in a Fabry-Perot cavity [51]. 
They suggested possible applications to solid-state physics 
and quantum field theory [52]. Remaud and Hernandez [53] 
found that a time-dependent mass parameter offers a means 
of simulating input or removal of energy from the system. 
Implementation of controls on these time-dependent quan- 
tum systems requires guidance from mathematical studies of 
controllability for time-dependent Hamiltonian operators. Al- 
though the HTC theorem deals with controllability in infinite- 
dimensional Hilbert space, it is restricted to time-independent 
operators. This paper explores a more general case. We seek an 
extension of the HTC theorem that is applicable both to time- 
independent and time-dependent quantum systems, as well as 
to systems with discrete or continuous operators acting on 
finite- or infinite-dimensional state spaces. 

Since this paper is aimed at an interdisciplinary readership 
that includes pure quantum theorists as well as control en- 
gineers, it is well to draw a clear distinction between time 
dependence of the system arising solely from influences that 
are directly under the control of an external, purposeful agent, 
and time dependence that is intrinsic to the physical system 
either in isolation or as embedded in a natural environment. In 
the accepted terminology of control theory, which we adopt, 
the former case defines a time-independent control system, and 
the latter, a time-dependent system. The issue of controllability 
has received considerable attention in the time-independent 
situation so identified (e.g., in Refs. [5], [8], [22], [12]); 
whereas relevant results for the time-dependent case are very 
limited. 

The time-dependent quantum control problem that we shall 
address is stated formally in Sec. 2. To cope with the unbound- 
edness of operators involved in the Schrodinger equation, an 
analytic domain is introduced in Sec. 3, such that solutions of 
the Schrodinger equation can be expressed globally in expo- 
nential form on this domain. In Sec. 4, we define an augmented 
system in a space enlarged by one dimension, enabling its 
description within the framework of time-independent control 
systems. Following the pattern of Kunita's proof [54] of 
strong controllability of a time-independent system, we then 
establish conditions for controllability of this kind for the 
one-dimension-reduced system defined by the original time- 
dependent Schrodinger equation. Three illustrative applica- 
tions of the theorem are presented in Sec. 5, and our findings 
are reviewed in Sec. 6. 



II. Problem Formulation 

The following quantum control system is derived by ap- 
plying the geometric quantization method [55] to a classical 
bilinear control system [56], [31]: 



H' (t) + (t) 



(2) 



^(*o) = "00 ■ 



Here, H' (t), and the H[if) with I = 1, 2, . . . , r, are Hermitian 
operators on a unit sphere Su of Hilbert space, the ui(t), I = 
1, . . . , r are restricted to piecewise-constant real functions of 
time, and ip(t) denotes a quantum state belonging to Su- 
In physical language, Hq is the unperturbed or autonomous 
Hamiltonian, and the H[ are interaction Hamiltonians. It is 
the coefficients ui(t) that are subject to purposeful control by 
an agent external to the system, within the specified class of 
functions. Setting h = 1 and dividing H' (t) and the H[{t) by 
i, we arrive at a more familiar control form, 



d_ 
dt 



(3) 



where the Hi(t), i = 0, 1,2, ...,r, are skew-Hermitian op- 
erators on Su- From the standpoint of systems engineering, 
Ho(t) is called the drift term in Eq. Q because no control 
function directly modifies its action. Importantly, we depart 
from previous studies of quantum controllability in allowing 
the Hamiltonian operators Hi(t) to their own carry explicit 
time dependence, which is assumed to be inherent in the 
physical structure of the system and therefore beyond the 
control of any external agent. The operators Hi(t) are the 
counterparts of the structural matrices involved in standard 
formulations of linear control theory. 

For the system Q, we know from arguments presented 
in Ref. [5] that the transitivity of states on Su requires an 
infinite sequence of control manipulations within the control 
set {ui(t)} of piecewise-constant real functions. Clearly, such 
a process is strictly meaningless in practice, although under 
certain conditions it may be possible to find a finite series of 
control operations that approach the desired target state arbi- 
trarily closely. Even so, we are naturally directed to consider 
the issue of controllability on a. finite-dimensional submanifold 
of the unit sphere Su, for which in turn a finite-dimensional 
tangent space is generated by H (t)ip(t), . . . , H r (t)ip{t). 

Accordingly, our attention focuses on a finite-dimensional 
submanifold M C Su, on which the following dynamics 
prevail 



d 



dt 



H (t) +J2 u l(t)H t {t) 



V> , Ht) € AT, Vfc > to : 



(4) 



Thus, instead of studying controllability on Su, we consider 
controllability on M C Su- On the submanifold M, the 



inherited topology of Su still applies; hence it is paracompact 
and connected. 

For system @, we have available a set of vector fields 
0(M) composed of skew-Hermitian operators on M with 
Lie algebra defined by O(M) = £{H , . . . , H r }. Let V be a 
subset of 0(M). The Lie algebra generated by V is denoted 
by C{V). The restriction of C(V) to a point ip on M, which 
is a tangent subspace of TM$ at ip, is written as 



while 



C{V){%I>) = {Y^\Y e C(V)} C TM^ , (5) 



C(V) = {£(V)V# e M} (6) 



defines an involutive differential system. A vector field X is 
said to belong to C(V) if X(ip) £ C(V)(i/j) holds for all 

ip e M. 

III. Selecting the Domain 

Recognizing that operators in quantum mechanics are in 
general unbounded, we need to find a domain on which 
exponentiations of the operators entering the system @ con- 
verge. To this end, we introduce the so-called analytic domain 
conceived by Nelson [32], a dense domain invariant under 
the action of the operators in system (|3j. The solution of the 
Schrodinger equation can be expressed globally in exponential 
form on this domain, which is also invariant under the action 
of the exponentiations of the operators Hi. 

Definition 3.1: If H is an operator on the state space 7i, 
we call an element uo of H an analytic vector for H in case the 
series expansion of cxp(Ht)uj has a positive radius of absolute 
convergence, that is, provided 



oc 

E 

n=0 



\H n uj\ 



-s < oo 



(7) 



for some s > 0. 

If H is a bounded operator, then every vector in Ti is trivially 
an analytic vector for H. 

The corresponding definition of analytic vectors for a Lie 
algebra of operators runs as follows [32], [57]: 

Definition 3.2: A vector u> e Ti is said to be an analytic 
vector for the whole Lie algebra C if for some s > and 
some linear basis {Hi, . . . , Hj} of the Lie algebra, the series 



1 y 



■ H t u\\s r ' 



(8) 



l<ii, 



is absolutely convergent. 

The concept of analytic vectors is especially useful for our 
purposes, since for certain types of unbounded operators they 
form a dense set in the Hilbert space. In fact, the set of all 
analytic vectors for a Lie algebra L forms an analytic domain 
in the following sense [32], [57]. 

Definition 3.3: Let L be the Lie algebra generated by the 
skew-Hermitian operators Hq, ■ ■ ■ , H r on a unit sphere Su of 
Hilbert space. An analytic domain T>a is said to exist for the 
H^ % = 0, 1, . . .r, if (i) there exists a common dense invariant 



subspace T>a C H on which the corresponding unitary Lie 
group G can be expressed locally in exponential form with 
Lie algebra £, (ii) T>a is invariant under G and C, and (iii) on 
T>a, elements of G can be extended globally to all f € M. + . 

We now state Nelson's fundamental theorem, which pro- 
vides conditions under which a Lie algebra C defined by a set 
of skew-Hermitian operators can be associated with a unitary 
group G having L as its Lie algebra. 

Theorem 3.1: (Nelson) Let £ be a Lie algebra of skew- 
Hermitian operators in a Hilbert space H which have a 
common invariant dense domain T>a- Let X\,.,. , be an 
operator basis for C. If T = X\ + . . . + X\ is essentially 
self-adjoint, then there is a unique unitary group G in H with 
Lie algebra C. Let T denote the unique self-adjoint extension 
of T. Then the analytic vectors of T are analytic vectors for 
the whole Lie algebra C and form a set invariant under G and 
dense in 7i. 

Accordingly, on the analytic domain T>a, the Lie algebra 
and its unitary Lie group are related through the familiar 
exponential formula. The Lie algebra is composed of skew- 
Hermitian operators which are vector fields defined on T>a H 
Sn- By property (iii) of the definition 13.31 of the analytic 
domain, these vector fields on T>a H Sn are complete. More- 
over, owing to the skew-Hermiticity of the operators Hi of 
system ©, the corresponding transformation groups, taking 
a given state on Sn to another state on Sn, are unitary. 
This feature guarantees preservation of the norm of quantum 
states, as required for the statistical interpretation of quantum 
mechanics. 

In fact, Nelson's theorem only provides sufficient conditions 
for the important properties it yields. With this in mind, we 
shall assume an analytic domain T>a exists without explicitly 
imposing the conditions stated in this theorem, a stance also 
adopted in Ref. [5] This strategy clearly implies that the 
existence of such a domain must be established explicitly prior 
to application of the controllability results to be derived in the 
following sections. 

We are now prepared to adapt the concept of controllability 
to problems involving unbounded operators. 

Definition 3.4: For system if T>a exists for £, and 
if for any i/jq and ipf € T>a H Sn there exist control 
functions U\(t), . . . , u r (t), and a time tf [resp. Vf/] such 
that the solution of control system satisfies ip(to) — tpo, 
tp(tf) = ipf, and ip(t) £ V A H Sn, where t < t < tf, then 
the system is called analytically controllable [resp. strongly 
analytically controllable] on Sn', moreover we then say that 
the corresponding unitary Lie group is analytically transitive 
on Sn- 

As has been argued, the more pertinent concept is con- 
trollability on the submanifold M of Sn- By assumption, 
MnV A is dense in M, while dim(Mnr>A) = dimM = m. 
Denoting the tangent space of M n T>a at ij) by TM$ = 
C{Hq, . . . , H r }ip, the tangent bundle of the system l|4} is given 
by T(M n V A ) = U^Mn^TM^. 

Let Rt(ip) denote the set of all points that are reachable 
from ijj at time t. The set R(i/j) = Ut>t ^t(ip) is then 



reachable from ip at some time greater than to. We say 
that system (|4) is analytically controllable on M if R(ip) — 
M n T>a, V-0 € M n V A , and that the system is strongly 
analytically controllable on M if R t (ip) = M n T> A ,\/t > 

t , w e M nv A . 

IV. Controllability of Time-dependent Quantum 
Control Systems 

A. Reformulation as a Time-independent Augmented System 

Most of the methods developed for determining controlla- 
bility of time-independent bilinear or nonlinear systems [58], 
[59], [5], [31], [60], [61] cannot be applied directly to the 
time-dependent bilinear control problem studied here, since 
these approaches rely upon the following property. Let Y t (ip) 
be an integral curve of the time-independent tangent vector 
Y starting from point and t 6 [to, to + tf], an d l et 
cY t (ip) be an integral curve of the tangent vector cY starting 
from tp and t S [to, to + t//||c||]; then the integral curves 
Y t ((p) and cY t (ip) coincide. This property holds for all time- 
independent tangent vectors, but it generally fails for time- 
dependent tangent vectors. 

However, recognizing that this feature has been instrumental 
to controllability proofs for nonlinear systems, we recast the 
system @ as a time-independent problem so that it can once 
again be exploited. Reformulation of the original problem is 
accomplished by regarding the time variable t as an additional 
parameter in the specification of the system state, supplement- 
ing the state vector if). Thus the state of the extended system 
is expressed as 

' t + t 



(9) 



Making the corresponding extension of the manifold M, we 
form an augmented (m+ 1) -dimensional manifold defined by 



N = 



M<~)T>a 



(10) 



where R is the real line. Next we define augmented vector 
fields Wi by 



w (0 



1 



(11) 



£Zo(t + toMt + t ) 


Hi(t + t )i>(t + t ) 

with I = 1, 2, . . . , r. Obviously, the Wu with I = 0, 1, . . . , r, 
depend on both t and ip, i.e., the Wi now depend on the state 
£ defined by Eq. (|9j. 

The time-dependent control system @ has thereby been 
reformulated as an augmented system of time-independent 
form. Explicitly, 



dt 

£(0) = r? = 



Wo(£)+X>(«£) 



(12) 



to \ = I t 

vt > o, Vo e M nx>A, £ e N, 



where N is the n = (TO+l)-dimensional manifold constructed 
in Eq. dlOi and M is now viewed as a one-dimension-reduced 
manifold of the augmented system. As always, the controls 
ui(t), with / = 1, . . . , r, are piecewise-constant real functions 
of time t. 

It is convenient to employ t + to instead of t in definitions 
(|9jl and il y , thereby setting the starting time at zero for 
the augmented system dl 2b . Since the latter system is time- 
independent by construction, this can be done without affect- 
ing its trajectory. Thus, if the time for the augmented system 
is t, then the time for the original system 10} is t+t . Standard 
differential equation techniques can evidently be employed to 
analyze the behavior of the augmented system on the manifold 
N, and the results will reflect the behavior of the original 
system on manifold M, 

We note peripherally that system H2i is in a decomposed 
form in the sense of Ref. [59], where several theorems were 
developed for decomposition of nonlinear control systems. 
However, these theorems do not specify reachable sets, so they 
cannot be applied here to obtain controllability results. 

Reachable sets Rtiv) an d R(v) m ' e defined for the aug- 
mented system i\2i in just the same manner as for system 
@. From the work of Huang, Tarn, and Clark [5] based 
on the results of Chow [62], Sussmann and Jurdjevic [24], 
and Kunita [54], [58], it is to be expected that the issue of 
analytic controllability will hinge on the relationships among 
certain Lie algebras generated by the vector fields involved 
in the control system @ or its augmented counterpart dl 2i . 
For the latter problem, these Lie algebras are specified by 
A = C{W ,...,W r }, B = C{Wt,...,W r }, and C = 
£{ad^ Wi, I = 1, . . . ,r, m = 0, . . . , oo}. By definition, 
ad^ Wj is built from repeated commutators of Wo, present 
in A but not B, with any and all of the Wi present in A or 
B; clearly, 

BcCcA. (13) 

For future reference we note (in particular) that the restriction 
of B to a point ip on TV, which is a tangent subspace of TN^ 
at ip, is written as 

B(ip)={Y(iP)\Y eB}cTN^, (14) 

and in turn that 

B = {B{ip)\ip G N} (15) 
is an involutive differential system. 

B. Controllability of the Augmented System 

We must still face the situation that standard controllability 
results [58], [59], [5], [31], [60], [61], derived for time- 
independent systems, cannot be carried over directly to our 
problem as reformulated in the preceding subsection, since 
derivation of these results employs the vector-space property 
of the tangent space. Specifically, it is required that if Y is 
an acceptable tangent vector, then so is cY, where c is an 
arbitrary constant. But in our case, once the first component 
of a tangent vector of the augmented manifold is fixed at 



unity, it is not possible for both Y and cY, with c ^ 1, to 
be available tangent vectors. However, with the aid of a result 
of Kunita [54], we may nevertheless establish one-dimension- 
reduced controllability of the augmented system; that is, we 
may prove strong analytic controllability of the original system 
since it is not necessary to control the time dimension. 

First, let us identify certain properties of the reachable 
set Rt(?l) that will be useful in proving strong analytic 
controllability. 

Theorem 4.1: [24], [54] Assume that the Lie algebra C 
is locally finitely generated, and let I(rj) be the maximal 
connected integral manifold of C containing the point rj. Then 
Rt(v) C att(I(r))), where a® is the integral curve whose 
vector field is Wo. Furthermore, the interior of Rt(i]) with 
respect to the topology of a^(I(rj)) is dense in Rt(v)- 

A key relationship between the interior of the reachable set 
Rtiv) °f the augmented system at time t and the interior of 
its closure is provided by the following lemma. 

Lemma 4.2: 

int (cl 4>) )= int £*(??). (16) 
Proof: Let x G int(cli? t (?/)) and let S e (x) be the set of 
all \' such that x is reachable from x' within time e > 0. 
Then S e (x) is the reachable set within time e > for the dual 
control system 



Or 



W (v) + J2ui(t)Wi(v) 



(17) 



Theorem 4.1 implies that intiS e (x) is dense in c\S e (x), and 
int Rt(rj) is dense in cl Rt(rf). Since x G clS e (x)> we know 
that 

clS e ( X )n int (cl /Jtfa)) ^0 (18) 



and hence that 



int S e ( X ) n mt(cl R t (v)) H R t (r)) ^ . 



(19) 



If ( belongs to the latter intersection, then ( is reachable from 
7] using time t, and x is reachable from £ in elapsed time less 
than or equal to e. Therefore, x is reachable from 77 in elapsed 
time between t and t + e. This argument holds for any t > 
and any e > 0. Letting e — > 0, we conclude that x is reachable 
from t] in time t, so x G Rt(v)- Thus, 

int (cl 40/0) C Rt{v) => int (cl R t (77)) c int 4 fa) ■ 

But clearly int Rt (77) C mt(cl Rt(r/)) and the statement d!6i 
follows. 

From the control-theoretic perspective, the drift term is un- 
desirable because no control is present to influence or remove 
its effect. It is therefore of strategic value to consider a suitably 
modified control system, called the auxiliary system, that will 
serve as a bridge to an effective controllability analysis of 
the augmented system. Let eo,ei,...,e r be unit vectors in 
E r+1 ; in particular, let e t = (0, . . . , 0, 1, 0, . . . , 0), in which 
only the (i + \) th element is unity and the others are zero. 
Denote by Uq the set of controls u(t) — (uo(t), . . . , u r (t)) 



composed of piecewise-constant functions itj (t) taking the val- 
ues eo,±ei, . . . ,±e r only. Consider then the control system 
expressed in the form 

||=uo(t)Wo(0+5>i(WO, C(*o)=»>. (20) 
i 

where u(t) S Wo- The solution of this system may be written 

as 

(21) 



1 k l 3 ll 

a t = a t k ■ ■ ■ a t 3 ■ ■ ■ a t \ , 



where k is a positive integer and where a t J . is the integral 



curve of W; with u 



0,1, 



■,r, J 



, , k, and k a 



positive integer. The times tj satisfy tj > if ij — 0, tj € K. 
We denote by (77) the reachable set of the auxiliary system 
corresponding to the total time t since time zero, over which 
the control function ito(-) is nonzero; the reachable set of the 
auxiliary system is then R ^) = [J t>0 Rtiv)- Theorem 14.11 
is valid for this control system [24]. 
The following notations are convenient: 

Exp £ — the group of diffeomorphisms generated by 
the a\, t € K, i = 0, . . . , r, where a] is an 
integral curve of Wi , 
(Exp£) + = the semigroup of diffeomorphisms generated 
by a°, t > 0, and the with (eK 
and Z = 1 , . . . , r , 

(Exp£)t = the subset of (Exp£) + generated by 

fc 

a tt ■ ■ ■ ■ ■ a \\ , with *J ' 1 (b=o} = 1 ■ 

i=i 

To clarify the meaning of the last line, we note that when the 
index j is such that ij = 0, we have uq = 1 (and all the other 
Uj = 0), so Wo is "turned on" and does play a role as an 
active vector field or tangent vector. Conversely, for indices 
j such that ij ^ 0, the factor u multiplying Wq in system 
j20t vanishes, and Wo plays no role. The sum appearing in 
the definition of (Exp£) t gives the total time over which Wq 
is active in the system dynamics. 

From Chow's theorem [62], [24], it is known that the 
group Exp£ acts transitively on the manifold TV when 
dim £{Wo,Wi, ... ,W r } = dimiV, i.e., we know that 
{a(rf)\a £ Exp£} = TV for any 77 G TV. On the other hand, 
the reachable set at time t for the auxiliary system \2Q\ is 
R^i]) = {a(ri)\a <E (Exp£) t }. (It is to be noted that in the 
present context t is the total time over which Wo has been 
active since time zero, which is generally not equal to the 
actual elapsed time, since Wo may be turned off over certain 
intervals.) 

Lemma 4.3: 



c\R t (f 1 )=c\R° t ( V ). 



(22) 



We may gain intuitive understanding of this lemma by ana- 
lyzing a simple example. 



(23) 



(24) 



Example. Let us compare the control system 

^(y) = (o) +u (i 

wherein iieR, with the system 

Jt{l ) =Wo (o) +Ml ( ! 

wherein (uo, ui) € {(0, ±1), (1, 0)}. Clearly, the first of these 
corresponds to the augmented system, and the second to the 
auxiliary system. Let Rt{r)) and Rt(v) denote respectively 
the reachable sets of systems fl23i and ( 1241 . staring from the 
state rj. While stopping short of rigorous argument, explicit 
computation will be used to reveal the pertinent relationship 
between c\R t (rj) and cli?((7?). 
First consider the integral curve 

and for n = 1, 2, 3, . . . form a series of integral curves /3™(f?) £ 
R t (77) defined by 



As 11 goes to 00, we find 


1 









-1 



(26) 
(27) 



that is, /3"(?7) — > ctt(v)- Hence ott(rj) £ cl Rt(rj). 
On the other hand, consider 



ft fa) 



m 2 





-1 



1 



(28) 



where mi, mi £ M and t = t\ + ti + t%, and construct 



• mi 



(29) 



again for n — 1,2,3,.... Applying the Baker-Campbell- 
Hausdorff formula, it straightforward to show that 




(30) 



Similarly, let 



m 2 





-1 



(31) 



and employ the Baker-Campbell-Hausdorff formula to obtain 



lim ctg 

n — >oo 




(32) 



Obviously 



• <*l G £°(r?) , 



(33) 



and we find that 
lim < ( 

1 

o 



<2 



1 



- ni2 



mi 


-1 



ti 



t 3 



Therefore (3 t (v) & clB$(rf). 

Now let us proceed with the proof of Lemma 14.31 showing 
first that clEj(ry) C cl R t (r)). Consider that 0^(77) G RS(r]) is 
expressible in the form of a\ k k ■ ■ ■ a 1 ^ (77), where t = Ey=i tj ' 
l{i -=0}' With the guidance of the example above, a sequence 
of controls u^(-) associated with the diffeomorphism of this 
form is constructed as follows. For an arbitrary positive integer 
n such that nt m > J2i-^o \^j\> wnere m is the last subscript 
j such that ij — 0, let 



(n) _ t _ 



Define real numbers s\ 

< s[ n \ by 



» 



(35) 

n 

, sj™\ ordered so that < < 



< 



» _ 



(n) 
S J>2 



|*i I if *i = 0, 
£l*i| if ii/0- 

if last j with = , 

if other j with iy = , 



-^i + l^l 



» 



if i i9 4 0. 



(36) 



Further, let 

»(r) 



n • sgn(tj)e i;i . 







if -a 

if < r < Sj 



' r < sPtcij ^ 0, 



fn) 



& i ,■ = , 



if r > s 



(n) 



(37) 

where ei, . . . , e r are unit vectors in W. The solution of 
the system d!2l > associated with the control vS n '{-) may be 
written 

^)=cf^& e ^^' (38) 

where /3p|* J is the integral curve of Wo if ij = 0, or the 
integral curve of Wq + n ■ sgn(r)Wi j if ij ^ 0, i.e., 

P";i 3 = (W )r 



(W 



if ij = , 
sgn(r)W i3 )M if ij^O. 



(39) 



We note that (Wo + n • sgn(r)W / i .)M and {±W + 
sgn(r)M / i j . )| T | describe the same integral curve on N, by 
virtue of the time-invariance property of system d!2i . Obvi- 
a\ 3 as n — > 00. On the other hand, 



ously, (3 ]t 'J 

= <. 

Thus, as n 



j ' 1 {»'i=0} 



-> 00 we obtain 



(40) 



(41) 



and hence at(f]) G c\R t (rj). Because 0:4(77) is an arbitrary 
element in Ri(rj), it follows that (77) C c\R t (rj), and since 
cl R t (rj) is closed, it follows in turn that cl 1^(77) C cl J2 t (77) . 

Next we show c\R t (r]) C clijj(?7). Consider /3(r)) G Rt(v) 
of the form of • ... • /^(t?), with = expv,j(W + 



(34) c}Wi 



CjW r ) and Cj = (cj, . . . , dj). Here, cj- is the 



control applied to Wi during time period Uj, so Cj is the 
control set applied to Wi,...W r during the corresponding 



time interval u,-, with Uj G 



and c' G 



For each 



, fc, take a? in the form 



a 



Ui , 1 



7/ . 7^ ■ 

•exp^-(cW r )exp^-W 



(42) 



n J n J n 

Invoking the Baker-Campbell-Hausdorff formula [63], we 
write 



lim a™ 

n — >oo J 



= lim 



(43) 



exp ^-{c)Wx) ■ ■ • exp ^(<5W P ) • exp ^VK 



lim exp [«j(Wo + c)^ 1 



<£W r ) 



E 

0<p,g<r 



= exp Uj{W + c]Wi + ... + c r jW r ) = P c u ] ■ 
Constructing a™ . . . aJJ G R® (77) we then obtain 



lim 



*fc • ■ • 



(44) 



(45) 



so that /3(7/) G cli?° (77). Since /? (77) is an arbitrary element of 
Rt{ri), we arrive at Rt(rj) C cl J2j (77) and hence cl R t (r]) C 
cl-R°(r7). We conclude that cl^r?) = cl R^(rj). 

The time i labeling these reachable sets is to be interpreted 
as the time interval over which the control operation repre- 
sented by Wo is in effect, or "turned on." In fact, Wo is 
necessarily always "on" in the augmented system, so the total 
time elapsing in the augmented system is the same as the time 
interval over which Wo is turned on; hence the reachable sets 
R t corresponding to these two times are identical. Of course, 
the same coincidence does not hold for the auxiliary system. 
However, this is immaterial, since the auxiliary system was 
only introduced to exploit the key relationship J22b . Further, 
we may observe that the reachable set i?" (77) of system i20i , 
with the control u(t) = («o(t), ■ ■ ■ , u r (t)) assuming values 



(eo,±ei, . . . , ±e r ), is the same as the corresponding set for 
which the control u(t) assumes the values eo, ±cei, . . . , ±ce r , 
with ceR+. 

Since we can take advantage of the result J22b in this 
manner, it is clearly preferable to study the properties of 
Rtiv)- The auxiliary system is easier to control, and the state 
at time t can be expressed as a composition of integral curves 
of Wi in the same style as Eq. (12 1 1 . To do so, let the set 
of subscripts j with ij = be written as {p, . . . ,q,s} in 
increasing order, of course with t p + . . . + t q + t s = t. Then 
we have 



a t 



/ it >s+n / is-i \ i i 

i a t h ■ ■ ■ a C+\) ' ( a t ■ a t s . 1 ■ a -tJ ■ «, • 

„o 



K 



l <J-2 \ 



+*„H ht P 



a t P -i ■ a -(t s +t,H — ht P ; 

all • a 



,+*,+.. .+t P " "ti " "-(t a +t 9 +...+t p ) 

= po (an ) • • • A) (afc; ) • ft, (a£:l ) • ft, («*;: 2 2 ) • • • 
ft s +*,K:i) • Pu+tMl-l) • • -ftK:i) • • • 

•ft«)-a?> (46) 

where /3t(7) = a® ■ 7 • This analysis stimulates us to 
define the following three sets of diffeomorphisms: 

Exp/3 = the group generated by a\, t E R , Z = 1, . . . , r, 
where a| is the integral curve of vector field Wi , 

*i = U^=i {ft, (7*) • • ■ ■ • ft! (7i)l lj G Exp B, 
< ijfe < . . . < *i = t} , 

G t = U£i {ft fe (7fc) • • ■ ■ • A 1 (7i)l 7, G ExpS, 
ramtj > 0, maxtj = <}. 



By construction, 



A°(T 7 ) = f t a?(» J ) 



(47) 



We observe that F t is a semi-group of diffeomorphisms 
included in the the group Gt, whose properties are established 
in the following lemma. 

Lemma 4.4: First, the set Gf is a group. Furthermore, 
if dim C (77) = n — 1 = m holds for all 77 G iV, then 
{a(rj)\a G Gt} = a° (I(a°_ t (il))) is true f° r ah" ^ where 
is the maximal connected integral manifold containing 
v E N, whose associated Lie algebra is C. 
Proof: For a\,a-z E G t , it is easily seen that a\ ■ a<i G G t . 
Writing a G Gt as a = (3 tk (7^) ■ . . .■ (3 tl (71)1 we also see that 
= fti(7f 1 ) ' • ■ ■ ' ftfc^ 1 )- Therefore G t is a group. 

Now, denote the set {a(rf)\a E Gt} by B t (rj). It is straight- 
forward to show that (i) B t (r}) = B t (£) if £ G B t (r]) and (ii) 
£4(77) D -B t (£) = if £ £ B t (77) [54]. We can demonstrate 
that (iii) 77 G int .84(77) under the topology of a°(/(a^ f (77))) 
as follows. By definition, Rtiv) is me reachable set for the 
system ( 1201 . By the same reasoning that leads to Eq. (1471 . 
we have R°(a°_ t (v)) C ^(77) because R°(a°_ t (r])) = F t ■ 
a® ■ a Q _ t (rf). Since R^(a°_ t (ri)) has a nonempty interior with 



respect to the topology of aj (I(a^_ t (77))) by Theorem l4.ll we 
see that B t (rj) contains a non-null open set U. Given \i G U, 
choose a E G t such that 01(77) = Since a is a continuous 
map, a~ l {U) is an open set containing 7;. 

In fact, a~ l (U) is included in B t (rj). We know that Gt 
is a group, so oT Y E Gt if a € Gt. Letting £ e a _1 (t/), 
we can find x G U, such that x = a(Q E U C -84(77) 
and also \ E B t ((). By properties (i) and (ii), we obtain 
X G B t (QnB t (ri) ? 0. Hence B t (C) = ^(77) and C G ^(77). 
Accordingly, a _1 (J7) C -84(77) and 77 S int -84(77) under the 
topology of a°(7(a° t (j7))). 

The properties (i)-(iii) imply that Bt(rj) is maximally con- 
nected and open under the topology of a°(/(a^t(f7)))- Thus 
we have B t (rf) = a®(I(a°_ t (??))) for all t > and 77 e 
7Y. In addition, it is seen that £^(77) = a®(I(a Q _ t (ri))) = 

a r r^Ti I ■ The proof of Lemma I4~4l is now complete. 
Based on Lemmas 14.31 and 14.41 we could conclude that 

cl Rt(a°_ t (rj)) = a° (I(a° t ( 7 ?))) if we cou l d establish that 
F t — Gt- The following proof takes a slightly different path. 

Let Exp B denote the group of diffeomorphisms generated by 
all one parameter groups of transformations with respect to 
vector fields belonging to B. The sets F t and Gt are defined 
in the same way as F t and Gt, i.e. via Eq. (17), but with 
Exp B entering in place of Exp B. 

Obviously, F t C F t and Gt C Gt hold. We shall now 
establish that F t — G t . 

Lemma 4.5: Let X be a complete vector field belonging to 

B, and let jt be the one-parameter group of transformations 
generated by X. Assume \B,C](rf) C B(rj) is satisfied for 
all 77. Then df3 s (^ft) is an isomorphism between B(rf) and 
^{Ps(lt)(v)) f° r each 77, and F t = Gt is true for all t > 0. 
Proof: Since /3 s (-f tl ) ■ ft (7*2) = ft(7ti+t 2 ) holds, we have 
d/3 s (7 tl +t2 ) = rfft (7ti ) ' dft (7*2 ) • Hence it is enough to prove 
the lemma's assertion for sufficiently small \t\. Let Y tyS = 

d(3 s (j t )Z, where Z e B. For each value of s, /3 s (7t) with 
f e i? is the one parameter group of transformations generated 
by da° s X, while 

BY 

— f- = -dp.h t )[da a X, Z] = dj3 s (7t)[Z, da° s X] . (48) 
ot 

Therefore [Z, da a s X] E B by assumption, because da®X 

belongs to C = {C(t?)|t7 G iV} [64], [65]. 

Now we fix a point 77 of N and a value of s E R. Let 
Z 1 , . . . , Z n provide a basis of 6 in an open neighborhood 
U of 77. Then there exist G°° functions /y on Z7 such that 
[Z^da^] = holds in U. Let e be a positive 

number such that /3 S (7^(77) € J7 for \t\ < e, noting that 
13 s{lt) is a continuous map of t and /3 s (7o)( 7 ?) = 7 ?- Then 
d/3 s ( 7 t)[Z%da°X] = E"=i/i^ft(7t)^' for |t| < e. Set 
F J *(*) = dp s {^ t )Z j . Then ^(t), with |t| < e, satisfies the 
linear differential equation 



dW(f) 
dt 



/ j J] 



/^ fe (t)i = i, 



(49) 



The solution V J (t) can be written as V J '(t) = 
Ysk=i9jk[t)V (0), where (gjk) is a regular matrix. Also, 
we have V*(0) £ B(rf) and y fc (0 G 0(A>(7*)) (77). The map 
d/3 a (7i) : B(r?) -> B(/3 s (j t ))(v) is bijective because (5^) is 
a regular matrix. Moreover, dAs(7t) retains the structure of 
the Lie bracket with respect to da°X. This establishes that 
dp silt) is an isomorphism between 6(77) and B(A(7t ))(??) 
for |t| < e. Since 7^ = /3 8 (a) • 74 • /3 a (a0 -1 (with s fixed) 
is a one-parameter group of transformations generated by 

d/3 s (a)X and d(3 s {a)X belongs to B, we know 7^ (with 

t E R) belongs to Exp B. But Exp B is generated by all such 
7t, so we arrive at the relationship 



A (a) (Exp B) A (a) -1 C Exp/3, for aeg. 
Let a be any element of G t , written as 



Ai(7i)> k > 0, max*/ = t. 



(50) 



(51) 



By induction we can prove that there exist . . . , 71 of Exp B 
and < Sfc < . . . < Si = t such that 

(3 tk ( lk ) /3 tl ( 7 i) = P Sk (lk) ■■■■■ /3 Sl (7i). (52) 

Here we only consider the case k = 2. If *2 < ti, there is 
no need for proof. Suppose ^ > ti, and set £3 = £2 — £1. 
Then we may write A2CY2) ■ Ai(7i) = Ai(A 3 (72) ' 7i)- B Y 
relationship J50t . there exists 71 of Exp£> such that A3 (72) • 

7i ' A 3 (72) _1 = 7i» i-e-, A 3 (72) -7i = 7i ' A3 (72)- This 
implies 



PtM'Ptr (ll)=Ph (A3 (72)- 7l) 

= A 1 (7i)"A 2 (72) 



Ai (71 -As (72)) 



(53) 



More detailed proofs may be found in Refs. [54], [66]. 

Theorem 4.6: Suppose that dim C(rf) = n — 1 = m holds 
for all r) E N, and suppose that [B,C](rj) C 6(77) holds for 
all 77. Let 1(77) be the maximally connected integral manifold 
containing 7/ whose corresponding Lie algebra is C. Then 
a?(Ifo)) = R t (ri). 

Proof: Clearly we have {aa^(r])\a E F t } C {aaj (r])\a E 
F t }. In fact, the closures of these two sets coincide. Since 
Ft — Gt D Gt, it is seen that 

clR^v) = c\{aa t {r,)\a E Fj} 
= cl{aa°(?7)|aG Ft) 
= cl {aa° t (77) I a E G t } (by Lemma IO 
= cl a? (J(a° t (a? (77)))) (by Lemma l4~4l 
= cla?(/(7,)). 

(54) 

But Lemma 1431 tells us that cl (77) = c\R t {rj), so we 
obtain c\Rt(i]) — cla® (IOt])). And from Lemma [TBI we know 
that int Rt{if) = int(cl Jit (77)), which implies inti? t (?i) = 
0^(7(77)) under the topology of a° (1(77)). Finally, i? t (77) C 
oij (7(77)) by Theorem l4.ll and we arrive at At (77) = a° (I(rj)). 



C. Strong Analytic Controllability of the Actual System 

In subsection 4.2, we investigated the reachable set at 
time £ of the time-independent augmented system formed 
by enlarging the state space to include an extra dimension 
corresponding to the variable t. Now we return to the original 
quantum control system @ to discover conditions under which 
it is strongly analytically controllable. 

Theorem 4.7: For the control system defined by Eq. 

let 

B(t)=C(H 1 (t),...,H r (t)) 



B X = -[H Q ,B] + ^B 



B n — — [Hq, B n -x] + -QtB n 



(55) 



C = £{/3,5 1; ..., £?„,...}. 

Suppose dim C(t)ij)(t) = m holds for all ip £ M n V A , 
and [B, C](t) c Bit) is the case for all t. Then the time- 
dependent quantum control system @ is strongly analytically 
controllable. 

Proof: We apply Theorem 14.61 to the augmented control 
system 

G3- To do so, we need to examine the Lie algebras 
B and C for this problem. For B we readily find 



B=£{Wl,...,W r } 



Hi(t) 







H r (t) 



£{Hi(t),...,H r (t)} 
Next let us construct C. For any 







W{rj) = W{t,ip) 
where 77 E N, we have 
&d Wo W= [W ,W] = 







Gfl, 



(56) 



(57) 



8 







d(t,1>) 




1 



S 





H{t)m 
1 

H (t)m 



Similarly, 



#o(£)V(*) 


-[H ,H]+dH/dt 







Setting i?i = — [i?o,0] + dB/dt, we may then derive 
ad^ o /3 = adH/ advKoB 

" W 

SiVW / V -[flo,-Bi]+0Bi/0t 



(58) 



(59) 



ad 



(60) 



Continuing in this fashion with 



Br, 



[H ,B n ^}+dB n ^/dt 



for n = 2, 3, . . ., we find 



ad Wo B 





B n tp(t) 



(61) 



(62) 



= £ 



Thus 

C = £{i3, ad Wo t3, ad iVo S ' ■ ■ •} 

° 1 ( ° ■ 

\ (63) 

£{B(t),B 1 (t),...,B n (t),...}4>(t) 


mm 

From the assumption that [B,C](i) C V(t), we have 
[B,C](^(i)cB(^(t),V(t). (64) 



Hence 











c 





Bil) 



(65) 



since o£(/(a° t (77))) 



, we obtain a^{I{rj)) 



so that [B,C]{rj) C S(r?), Vr; € AT. 

By assumption, dim C(t)ip(t) = m, Vip € H £>a, which 
implies that dim C(rf) = m = n — 1 holds for all rj € N, 
According to Theorem Pk6l a°(I(r])) = Rt(v), Vt > 0, and 

to 

v m n 2?a 

i + t 
m nv A 

Let 7r : N — > M fl 2?^ be the projection map that in effect 
annihilates the time-dimension of the augmented problem 
corresponding to the variable i, and brings us back to the 
original control system. In fact, the extension and projection 
maps mediate a one-to-one correspondence between the states 
of the augmented system and those of the original system. The 
simplicity of this relationship stems from the fact that t is a 
strictly increasing variable. 

To reiterate our strategy: We have dealt with the explicit 
time-dependence of the original control problem by adding 
an extra dimension to its state space, such that, as viewed 
in the augmented space, the augmented control problem is 
time-independent. After analyzing controllability within this 
extension, the results are projected to the original space by 
removing the extra time dimension, recovering the exact states 
of the original system. 

Accordingly, 7r(a? (I(rj))) = M n V A , while 7rR t (r]) = 
Rt+tjip), V-0 € M n V A . Hence R t (*ip) =Mf] V A ,Mt > t , 
and the system (0} is strongly analytically controllable on M. 

We may note that upon introducing the Lie algebra A(t) = 
C{Ho(t) , H\(t) , ■ ■ ■ , H r (t)}, it is readily established from 
property <[l3j that B C C C A for all t. 



To complete the formal analysis, we state two corollaries 
that devolve immediately from Theorem I4.7I 

Corollary 4.8: From the operators Hi entering control 
system @, form the Lie algebras B = £{Hi,...,H r } 
and C = £{B, a,dff B, . . . , ad^ o £>, . . .}. Suppose that the 
Hi do not possess explicit dependence on the time t, that 
dim Cip(t) = m holds for all ip e MfXD A , and that [B, C] C B 
is satisfied. Then the time-invariant system @ is strongly 
analytically controllable. 

Corollary 4.9: For the control system (0}, form the Lie 
algebra Bit) = £(Hi(t), . . . , H r (t)), and suppose that 
dim B(t)ip(t) — m holds for all i/j S M n T> A . Then system 
|4]is strongly analytically controllable. 

The latter corollary follows because [B,C](t) C B(t) must 
hold, once dim B(t)i/j(t) = m. 

V. Examples of Strong Analytic Controllability 

In this section, we present three examples that meet the 
criteria for analytic controllability enunciated in Theorem 4.7. 
The examples selected are relevant to problems of interest in 
mathematical physics or engineering applications of quantum 
mechanics. 

Example 1 The strong analytic controllability theorem can 
be applied to the simple degenerate parametric oscillator, a 
problem of importance in physics and engineering. Introducing 
an appropriate effective Hamiltonian allows the corresponding 
control system to be written in the form [67] 



= L{t)cfa + l - X {t) [e- 2<w V) 2 + 



(66) 

Here and a represent, in turn, the creation and annihilation 
operators of the pump mode of frequency uj(t), while 
is the time-dependent coupling function related to the second- 
order nonlinear susceptibility of the pumped medium. We may 
consider u>(t) and x(t) as control functions playing the role 
of the ui in Eq. (|4j, since they are real and can be adjusted 
to piecewise-constant functions of time t, outside the system 
itself. 

Following precedent [68], [69], [70], [71], we define the 
operators 

A+ = i(at)2 j X_ = i a 2 , K = i( a t a + aa t) ) (67 ) 

which satisfy the commutation relations of SU(1, 1), thus 

[K , K±] - ±K ± , [K+,K-] = -2K . (68) 



Setting 



H 
Hi 

Ho 



-iK , 



t K+ + e 2iut K_] 



[ e - Jwt K+ - e Iwt KJ\/2 



(69) 
(70) 

(71) 



the control system J66l > may be written in the more familiar 
form 

-V=[w(t)ffo+x(*)ffi(*)]V'- (72) 



The skew-Hermitian operators Ho, Hi, and Hi satisfy the 
commutation relations 



[Hq, Hi 



-Hi , \Hq,Hi\ — Hi , [Hi,H 2 ] = H . 



(73) 

We observe that the system i72i does not have a drift term 
in the usual sense, because the factor oj(t) can be manipulated 
externally. We also see immediately that A = B = C = 
C{Hq, Hi, H2}, and the second condition of Theorem 4.7 is 
obviated. In addition, Hq has eigenvectors \mk), with rn = 
0, 1, . . . and k — 1/4, 3/4, which span an analytic domain 
T>a [69], [71]. Consequently, we can choose a manifold M 
such that dim Ci/j = dim M \/tp £ T>a f~l M. All conditions of 
Theorem 4.7 being met, the system (I66> is strongly analytically 
controllable on M. 

Example 2 Defining Q = idt + d XlXl + d X2X2 , the 
Schrodinger equation for a free particle moving in two spatial 
dimensions may be expressed simply as Qu = 0. Determina- 
tion of the maximal symmetry algebra of this equation leads 
to the following set of nine operators, which form the basis 
of a nine-dimensional complex Lie algebra: [72] 

K 2 = -t 2 d t - t{xid Xl + x 2 X2 ) - t + (t/4)(ar? + x%) , 
K-2 — dt, J — xid X2 - x 2 d Xl , Bj = -td Xj + ixj/2, 
Pj = d Xj , E = i, D = xid Xl + x 2 d X2 + 2td t + 1, 

(74) 

with j = 1,2. Of immediate concern is the real Lie algebra 
spanned by this basis, i.e., the Schrodinger algebra, which has, 
as alternative basis, the operators Bj, Pj, and E (yielding the 
five-dimensional Weyl algebra), plus the operator J and the 
three operators defined by Li = D, L 2 = K 2 + K- 2 , and 
L3 = K- 2 — K 2 . The pertinent nonvanishing commutators 
are specified by [72]: 

[Li,L 2 ]=-2L 3 , [L 3 ,Li]=2L 2 , [L 2 ,L 3 ]=2L 1 , 
[Li,Bj] = Bj , [Li,Pj] = -Pj , [P j: J] = (-l) J ' +1 Pi , 
[Bj, J] = (-iy +1 B l , [L 2 ,B 3 ] = -Pj , [L 3 , Bj] = -P 3 , 
[L 2 , Pj] = Bj , [L a , Pj] = -Bj, [Pj,Bj] = E/2 , 

(75) 

where j,l = 1,2, j ^ I. 

Now we consider the controllability of the system 



_5 

di 



t/j= [L 2 + ui(t)Li+u 2 (t)L 3 + u 3 {t)Pi+u 4 (t)J]i;. (76) 



In this case there is a time-dependent drift term in the 
vector field driving tp. The relations (175 \ imply the equalities 
B = C = C{Li,L 2 ,L 3 ,Pi,P 2 ,Bi,B 2 ,J,E], while the 
required analytic domain T>a is furnished by the span of 
the eigenfunctions ip n ,m of L 3 . These take the explicit, time- 



dependent form [72] 

VVm = {2 m+n+1 TTn\m\)-^ 2 exp[iTr(m + n- l)/2] 

(m+n)/2 



x exp 



(«f iv 3 ) 



V3 + \ 
v 3 - i 



H m {vi/V2)H n (v 2 /y/2) 



v 3 



(77) 



where x\ = Vi{\ + v^) 1 ^ 2 , x 2 — v 2 {l + V 2 ) 1 / 2 , and t = v 3 . It 
follows as before that the system J76b is strongly analytically 
controllable. 

Example 3 A quantum control system with position- 
dependent effective mass m — (2Ax)^ 1 has been described 
by the time-dependent Schrodinger equation [73] 

. d 



[iBI +u 1 (t)A(t)I I-+iu 2 (t)C]ip, 



(78) 



where B, C € M and A(t) is a real function of time t but 
in general not piecewise-constant. The operators Lj and I±, 
which are independent of time, provide a basis for an su(l, 1) 
algebra, and have the concrete realization 



J_ = -d x , I =xd x + 1, 1+ = x 2 d x 4 
which satisfies the commutative relations 

[I ,I±] = ±I±, [I-,I+] =-2J . 



2x. 



(79) 



(80) 



This effective-mass problem arises in the study of semiconduc- 
tor heterostructures and, more generally, of inhomogeneous 
crystals [74]. In the semiconductor application, the effective 
mass of a carrier depends spatially on the graded composition 
of the semiconductor alloys used in the barrier and well 
regions of the microstructures [75]. 

The wave functions of the stationary states of Eq. d78l > can 
be written as 

ip E {t,x) = -^=expl-iE [ B(a)da+ [ [-C(cr) 
\2-k L Jo Jo 

-^B(a)]da^ x exp {-01 (t) (xd xx + d x )} x~ iE - 1 ' 2 

= --L exp j-iE J B(a)da + J [-C{a) - ^B{a)]da^ 

00 n 1 

J2U(iB(t)E + -+l) 2 [-ai(t)r 



„-iE-n-l/2 



n=0 1=0 



(81) 



These eigenfunctions span the analytic domain relevant to 
Theorem 4.7. 
Let us define 

H Q = Bio + u 2 {i)C, , Hi = -iA{t)IoI-, (82) 

where we take u 2 (t) = —B/2C. Eq. ( I78> can be recast as the 
control system 

d 



dt 



V> = [Ho + ui{t)Hi]il) . 



(83) 



Here the drift term is time-independent. Using the commuta- 
tion relations J80> . we obtain [Ho, Hi] = —BH±, Obviously, 
B = C C A, so [B, C] = B. Choosing a manifold M such that 
dimAf = dimCV> for all ip £ M, we are assured that system 
J78> is strongly analytically controllable. 

VI. Conclusions 

In this paper, we have formulated the time-dependent quan- 
tum control problem and studied its controllability. Acknowl- 
edging the unbounded nature of operators commonly involved 
in quantum control systems, our analysis has been predicated 
on the existence of an analytic domain [32] on which exponen- 
tiations of such operators are guaranteed to converge. Within 
this framework, we have extended the established treatment of 
time-independent quantum control problems by introducing an 
augmented system described in a state space that is enlarged by 
one dimension, yet embodies the true dynamics of the original 
system. With the aid of techniques and results developed by 
Kunita [54], [58], we are able to explicate the one-dimension- 
reduced controllability of the augmented system. Projection 
onto the original state space then yields a proof of the analytic 
controllability of the original time-dependent quantum control 
system, under conditions similar to those required in the 
time-independent case. The theorem so established has been 
illustrated with examples drawn from mathematical physics 
and systems engineering. 
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